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Abstract. The kinematical formalism for describing spinning particles developped 
by the author is based upon the idea that an elementary particle is a physical system 
with no excited states. It can be annihilated by the interaction with its antiparticle 
but, if not destroyed, its internal structure can never be modified. All possible states 
of the particle are just kinematical modifications of any one of them. The kinematical 
state space of the variational formalism of an elementary particle is necessarily a 
homogeneous space of the kinematical group of spacetime symmetries. By assuming 
Poincare invariance we have already described a model of a classical spinning particle 
which satisfies Dirac's equation when quantized. We have recently shown that the 
spacetime symmetry group of this Dirac particle is larger than the Poincare group. 
It also contains spacetime dilations and local rotations. In this work we obtain an 
interaction Lagrangian for two Dirac particles, which is invariant under this enlarged 
spacetime group. It describes a short- and long-range interaction such that when 
averaged, to supress the spin content of the particles, describes the instantaneous 
Coulomb interaction between them. As an application, we analyse the interaction 
between two spinning particles, and show that it is possible the existence of metastable 
bound states for two particles of the same charge, when the spins are parallel and 
provided some initial conditions are fulfilled. The possibility of formation of bound 
pairs is due to the zitterbewegung spin structure of the particles because when the 
spin is neglected, the bound states vanish. 
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1. Introduction and scope 

We introduce in the next section a brief description of the kinematical formalism for 
describing classical elementary spinning particles. The kinematical group of spacetime 
symmetries not only states the symmetries and conservation laws but also provides 
the classical variables we need to describe an elementary particle. The Dirac particle, 
a classical system which satisfies Dirac's equation when quantized, is summarized in 
section [3l It is Poincare invariant but the spacetime symmetry group has been recently 
enlarged [2] to include spacetime dilations and local rotations, so that the classical 
Dirac particle has a larger kinematical group of symmetries. Classical electromagnetism 
usually analyses interactions between spinless test point particles. If electromagnetism 
is a long range and also a short-range interaction it suggests that the same law applies 
at large and short distances and therefore that a scale invariance is involved. We have 
models of spinning particles which have a spacetime dilation invariance and we take 
advantage of this fact to analyse a long and a short range interaction among them 
and to check its behaviour when the spin is supressed. In section IH we construct a 
Lagrangian for describing the interaction between two Dirac particles which is invariant 
under this enlarged group. Because it is also invariant under spacetime dilations, it will 
describe a long and short-range interaction, which suggests a kind of generalization of 
the instantaneous electromagnetic interaction, because it has a Coulomb limit when the 
spin of the particles is smeared out. We also discuss there its structure in a synchronous 
time description of the evolution. Section [5] is devoted to the analysis of the evolution 
of two particles of the same charge when the spins are parallel to each other. One of 
the salient features is the possibility of formation, from the classical point of view, of 
metastable bound states for Dirac particles of the same charge, when some boundary 
conditions are fulfilled. 

2. Summary of the kinematical formalism 

Because all known elementary particles, the quarks and leptons, are spinning particles 
and it seems that there are no spinless elementary particles in nature, we took the 
challenge of obtaining a classical formalism for describing spin. The interest of a classical 
description of spinning matter is not important in itself, because matter, at this level 
behaves according to the laws of quantum mechanics. But finer a classical description 
of elementary matter a deeper quantum mechanical formalism, because we will have at 
hand, when quantizing the system, more classical variables to deal with, and therefore 
with a more clear physical and/or geometrical interpretation. A second feature is that 
a classical formalism supplies models. Both goals, in my opinion, have been successfully 
achieved. 

The kinematical formalism for describing elementary spinning particles [1], 
previously aimed for the spin description of matter, has proven to be a general framework 
for the description of elementary particles, because it supplies a very precise definition 
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of an elementary particle which has, as a quantum counterpart, Wigner's definition. 
All elementary systems described within this formalism have the feature that when 
quantized their Hilbert space of pure states carries a projective unitary irreducible 
representation of the kinematical group. It is through Feynman's path integral approach 
that both formalisms complement each other. 

It is based upon the four fundamental principles: Relativity principle, atomistic 
principle, variational principle and uncertainty principle. The relativity principle states 
that there exists a set of equivalent observers, historically called inertial observers, for 
whom the laws of physics must be the same. They are defined with respect to each 
other by a spacetime transformation group, usually called the kinematical group of the 
formalism. The atomistic principle admits that matter cannot be divided indefinitely. 
Matter does not statisfy the hypothesis of the continuum. After a finite number of 
steps in the division of matter we reach an ultimate object, an elementary particle. 
The distinction between an elementary particle and any other finite system is that an 
elementary particle has no excited states and, if not destroyed, it can never be modified, 
so that all possible states are only kinematical modifications of any one of them. This 
imphes that the states of an elementary can be described by a finite set of variables. The 
variational principle recognizes that the the action of the evolution of any mechanical 
system between some initial and final states must be stationary. This completes the 
classical framework. For the quantum description we must susbtitutc this last variational 
principle by the uncertainty principle, in the form proposed by Fcynman: all paths of 
the evolution of any mechanical system between some initial and final states arc equally 
probable. For each path a probability amplitude is defined, which is a complex number 
of the same magnitude but whose phase is the action of the system between the end 
points along the corresponding path. In this way, classical and quantum mechanics are 
described in terms of exactly the same set of classical variables. 

This formalism determines that these variables, which define the initial and final 
states of the evolution in the variational description, are a finite set of variables which 
necessarily span a homogeneous space of the kinematical group. We call them the 
kinematical variables of the particle. The manifold they span is larger than the 
configuration space and in addition to the time and the independent degrees of freedom 
it also includes the derivatives of the independent degrees of freedom up to one order 
less the highest order they have in the Lagrangian. The Lagrangian for describing these 
systems will be thus dependent on these kinematical variables and their next order 
time derivative. If the evolution is described in terms of some group invariant evolution 
parameter r, then, when writting the Lagrangian not in terms of the independent degrees 
of freedom but as a function of the kinematical variables and their r— derivatives, it 
becomes a homogeneous function of first degree of the r— derivatives of all kinematical 
variables. 

The formalism is completely general in the sense that it can accomodate any 
kinematical group we consider as the spacetime symmetry group of the theory. But 
at the same time it is very restrictive, because once this kinematical group is fixed 
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the kind of classical variables which define the initial and final states of an elementary 
particle in a variational approach, are restricted to belong to homogeneous spaces of the 
group. This kinematical group is the fundamental object of the formalism and must be 
defined as a preliminary statement. We call the formalism kinematical, to stress this 
fact. For the Galilei and Poincare groups, a general spinning elementary particle is just 
a localized and orientable mechanical system. By localized we mean that to analyse its 
evolution in space we have just to describe the evolution of a single point r, where the 
charge is located and in terms of which the possible interactions are determined. This 
point r also represents the centre-of-mass of the system for spinless particles, while for 
spinning ones must necessarily be a different point than q, the centre-of-mass, very well 
defined classically and where we can locate the mass of the particle. It is the motion of 
the charge around the centre of mass which gives rise to a classical interpretation of the 
zitterbewegung and also to the dipole structure of the particle. By orientable we mean 
that in addition to the description of the evolution of the point charge we also need to 
describe the change of orientation of the system by analyzing the evolution of a local 
comoving and rotating frame attached to that point. An elementary spinning particle 
is thus described as we use to describe a rigid body but with some differences: we have 
not to talk about size or shape and the point does not represent the centre-of-mass but 
rather the centre-of-charge. It is allowed to satisfy a fourth order differential equation 
and, for a Dirac particle, it moves at the speed of light. 

3. A Dirac particle 

For a relativistic particle the spacetime symmetry group is the Poincare group. In 
a recent work [2] we have shown that this spacetime symmetry group for a classical 
spinning particle which satisfies Dirac's equation when quantized, can be enlarged to 
include also spacetime dilations and local rotations of the body frame associated to 
the particle. We denote this group for the classical Dirac particle by iS = W ® S0{3)l, 
where W is the Weyl group, i.e., the Poincare group V enlarged with spacetime dilations 
and 5*0(3)^ is the group of local rotations of the body frame, which commutes with W. 
The Lagrangian for a free Dirac particle is also invariant under this enlarged group S. 

If we consider this new group as the kinematical group of the theory, then the 
kinematical variables of a Dirac particle are reduced to time t, position of a point r, 
where the charge of the particle is located, its velocity u with the constraint u = c, the 
orientation ex. which can be interpreted as the orientation of a local frame with origin 
at point r and characterized by three parameters of a suitable parameterization of the 
rotation group and, finally, a phase /3 of the internal motion of the charge around the 
centre-of-mass. If the particle has spin 5" 7^ and mass m 7^ 0, then a length scale factor 
R = S/mc and a time scale factor T = S/mc^ can be defined, such that all kinematical 
variables for the variational description can be taken dimensionless. It is this argument 
which justifies the enlargement of the spacetime symmetry group, to include spacetime 
dilations which preserve the speed of light. 
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The spinning particle has thus a centre-of-mass q, which is always different than 
the point r, such that, for a free particle and for the centre-of-mass observer is at rest, 
and the point r is moving in circles at the speed of light, around the point q in a flat 
trajectory contained in a plane orthogonal to the spin S. The kinematical variable (3 
describes the phase of this internal or zitterbewegung motion. 

The rotation subgroup of V transforms the kinematical variables r, u and ct among 
the inertial observers, and the local rotation subgroup S0{3)l only affects to the local 
change of the particle frame, i.e., to the orientation variables ex. Because the Weyl 
group has no central extensions [3], in the quantum representation the symmetry group 
becomes S = W® SU{2)l. 

The general structure of the Lagrangian for a free Dirac particle when written in 
terms of the kinematical variables and their r— derivatives is 

L = Ti + R-r + U-u + W-uj + B$ 

where T = dL/di, R = dL/dr, U = dL/dii, W = dL/du} and B = dL/dp, because 
of the homogeneity of L in terms of the r— derivatives of the kinematical variables. 




Figure 1. Motion of the charge of a Dirac particle in the centre-of-mass frame. 
The total spin S is the sum of the orbital part Z and the rotational part of the 
body frame W. It is not depicted the local body frame, with origin at point 
r, Cj, i = 1,2,3, which rotates with angular velocity cj. The motion of the 
charge, with respect to the fixed spin direction, is left-handed. The phase /3 
corresponds to the phase of the circular motion. 

The spin is defined by 

S = ux U + W = Z + W 

where Z is the zitterbewegung part of spin and W is the rotational part as depicted in 
Fig.lll The Z part is related to the function U, i.e., to the dependence of the Lagrangian 
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on the acceleration ii. When quantizing the system the Z part of the spin quantizes 
with integer values while the spin 1/2 is coming from the rotation part W. It is this 
twofold structure of the spin which produces a clear interpretation of the gyromagnetic 
ratio because the particle current, which produces the magnetic moment, is only related 
to the zitterbewegung part of the spin Z, which is twice the total spin S, as seen in 

The Casimir operators of the enlarged group S are the absolute value of the spin 
S, which is the Casimir operator of the Weyl group W and the absolute value / of the 
spin projection operator on the body frame of the rotational part of the spin 

which corresponds to the Casimir operator of the S0{3)l part. Here e^, i = 1,2,3 
represent the three unit vectors of the local frame attached to the point r. 

A Dirac particle, with the enlarged group S as its kinematical group, has as intrinsic 
properties the spin S and the spin projection / which take both the eigenvalue 1/2 
when quantized [2] . When quantizing the system only two kinds of Dirac particles arise 
according to the eigenvalues of the Z part of the spin, which can take the values or 1. 
The mass of the particle, very well defined as usual by means of the Casimir operator 
of the Poincare group — P^, and which must be different from zero in order to 
properly define the Casimir operator 5*^, is no longer an intrinsic property of a Dirac 
particle, because this operator does not commute with the generator D of space time 
dilations. Mass, at this stage with a larger spacetime symmetry group, represents a 
property associated with the different states of the Dirac particle. It could mean, for 
instance, that the electron and a massive neutrino are two different states of the same 
particle, the lepton. 

It is shown that the dynamical equation of point r for the free particle and in the 
centre-of-mass frame is given by 

r = -^S X u, (1) 

and where the spin vector S is constant in this frame, as depicted in Fig. [H The 
average value of the position < r >= qr is the centre-of-mass position, which is zero in 
this frame. The average value of the velocity < u >= v represents the centre-of-mass 
velocity, which is also zero in this frame. The radius of the zitterbewegung motion is 
R = S/mc, which becomes h/2mc, when quantized, i.e., half Compton's wavelength. 
This dynamical equation ([1]) is independent of the particular free Lagrangian we take 
provided the kinematical variables remain the same as the ones quoted here. 

4. An interaction Lagrangian 

An elementary particle can be annihilated by the interaction with the corresponding 
antiparticle, but if it is not destroyed, we made the assumption that the structure of an 
elementary particle is not modified by any interaction so that its intrinsic properties, the 
spin S and the spin projection on the body frame / cannot be altered by the interaction 
with an external field or by the presence in its neigbourhood of any other particle. 
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Let us consider a compound system formed by two spinning particles with the 
same kind of kinematical variables. We shall use a subscript a = 1, 2 to distinguish the 
variables corresponding to each particle. Then the kinematical space of the compound 
system is spanned by the variables {ta,ra,Ua,ota, Pa),o, = 1,2. The Lagrangian of the 
system will be written as 

L ^ Li + L2 + Li 

where the La, a — 1,2, arc the free Lagrangians of each particle and L/ is the interaction 
Lagrangian we are looking for. Both La are invariant under the enlarged group S 
and we are going to find an interaction Lagrangian Lj also invariant under S. The 
general structure of the free Lagrangian La of each particle, which only depends on the 
corresponding kinematical variables of particle a, is 

La = Tja + Ra-ra + Ua-Ua+Wa-UJa + Bja 

where Ta = dLa/dia, Ra = dLa/dVa, Ua = dLa/diia, Wa = dLa/dijJa and Ba = 
dLa/dj3a, because of the homogeneity of each La in terms of the r— derivatives of the 
corresponding kinematical variables. The spin and the spin projection on the body 
frame for each particle, are 

5„ = X ?7„ + Wa, hi = eai ■ Wa 

where Caj, i = 1, 2, 3 are three orthogonal unit vectors with origin at point Va- 

The interaction Lagrangian between these two particles Lj will be in general a 
function of the kinematical variables of both particles and of their r— derivatives. If 
both intrinsic properties Sa and /„ of each particle are not modified by any interaction 
then the interaction Lagrangian cannot be a function of the derivatives of the kinematical 
variables iia and uJa, a = 1,2. Otherwise the functions Ua and Wa will be different than 
in the free case. In this case the functions [/„ and Wa, which give rise to the definition 
of the spin, are obtained only from the corresponding free Lagrangian La- 

Then, as far as the r— derivatives of the kinematical variables are concerned, the 
interaction Lagrangian Lj will only depend on the variables ia, Va and a = 1, 2. In 
addition to this, it will also be a function of the kinematical variables ta, Va, Ua and 
j3a, but not of cx.a because of the invariance under the local rotation group 5*0(3)/,. 
Spacetime dilation invariance implies that the Lagrangian is a function of the phase 
difference (3i — P2, and of $1 — (52, but being both phases completely arbitrary and 
independent of each other it means that must be independent of these variables. 

Because of the constraint Ua — Ta/ia-i the interaction Lagrangian will thus be finally 
a function 

Li = Liita, Va, ta, To), 

and a homogeneous function of first degree of the derivatives ta, Ta, a = 1,2. 

If we call as usual the Minkowski four- vector x'^ = (ta,ra), translation invariance 
implies that the Lagrangian must be a function of — Xg. The following two terms 
Vfj.ui^i^2 r)^^{xi —X2){x2 — Xi), where rj/^i, is Minkowski's metric tensor, are Poincare 
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invariant. If we consider that the evolution parameter r is dimensionless, these terms 
have both dimensions of length squared. It therefore implies that its quotient is 
dimensionless and therefore invariant under spacetime dilations. Another requirement 
is that the Lagrangian is necesarilly a homogeneous function of first degree in the 
r— derivatives of the kinematical variables. The squared root will do the job. An 
additional discrete symmetry will be assumed because when the two particles are the 
same, and therefore indistinguishable, the interaction Lagrangian must be invariant 
under the interchange 1-^2 between the labels of the two particles. We thus arrive at 
the S group invariant Lagrangian 



' r^^,{xt - x^){xl - x\) (r2 - v^f - {t, - t,Y ^ > 

where gi is a coupling constant. This Lagrangian describes an interaction which is scale 
invariant and thus it is valid as a long and short range interaction and which has a 
Coulomb-like behaviour when the spin is supressed, as we shall see in the next section. 
In this way it suplies a kind of generalization of an action at a distance electromagnetic 
interaction between spinning particles. 

Alternative Lagrangians which fulfil the requirements of invariance and 
homogeneity can be constructed. For instance, L = g{xi — X2)'^{xi — X2)^/{xi — X2)'^-, 
but this one is a total r— derivative of the function log(a;i — X2Y . Another could be 
L = g{xi + X2)'^{xi — X2)^/ {xi — X2Y-, also dimensionless and linear in the derivatives of 
the kinematical variables, but it reverses its sign under the interchange 1 <-> 2, and thus 
all interaction observables, like the interaction energy are reversed, which is physically 
meaningless for two alike particles. 

The interaction Lagrangians L — gXiX2ix/{ii — X2)^{xi — ^2)^ and L — g\{xi + 
3^2)^(3^1 — X2)iu\/{xi — X2Y also fulfil the requirements of being invariant under the 
enlarged group <S, homogeneous of first degree in the derivatives of the kinematical 
variables, and invariant under the interchange of the particle labels, but they do not 
have a Coulomb-like limit when the spin is smeared out. They will not be considered 
here and the analysis of the interaction they describe is left to a subsequent research. 
Perhaps many other Lagrangians could be guessed according to the above requirements 
but we have not succeded in finding another one with such a desirable Coulomb-like 
behaviour. 

The interaction between two Dirac particles is not unique. We know that among 
leptons and quarks there are at least, short range interactions like the weak and strong 
interactions and a short and long range one like the electromagnetic interaction. The 
novelty is that this classical interaction Lagrangian, which possesses a Coulomb limit, 
describes the interaction between two spinning particles, which satisfy Dirac's equation 
when quantized, and has been obtained by assuming a spacetime symmetry group larger 
than the Poincare group. 



Interaction Lagmngian for Dirac particles 



9 



4-1- Synchronous description 

Once an inertial observer is fixed we shall consider a synchronous time description, i.e. 
to use as evolution parameter the own observer's time t which is the same as the two 
time variables ti and t2- In this case, t = ti = t2, ti = t2 = 1, and thus 




2 = 9 : (3) 



[r2 - ri)^ r 

where r = |ri — is the instantaneous separation between the corresponding charges 
in this frame and Ua = dva/dt the velocity of the charge of particle a. 

An average over the charge position and velocity in the centre-of-mass of particle 
1 implies that < ri >= qi and < ui >= 0, so that the interaction becomes the 
instantaneous Coulomb interaction, between the centre-of-mass of the first particle 
and the charge position of the other. The average over the other then corresponds 
to the instantaneous Coulomb interaction of two spinless point particles because when 
neglecting the zitterbewegung we are suppressing the spin structure. 

It is suggesting that g ~ ±e^ in terms of the electric charge of each particle and 
where the ± sign depends on the kind of particles either of opposite or equal charge. 

5. Analysis of a two-particle system 

The dynamical equation of a free Dirac particle and for any inertial observer is a fourth- 
order differential equation for the position of the charge r which can be separated into 
a system of coupled second order differential equations for the centre-of-mass q and 
centre-of-charge r in the form: [5] 



0, 



[q - r)2 



where now the dot means time derivative. The first equation represents the free motion 
of the centre-of-mass and the second a kind of relativistic harmonic oscillation of point 
r around point q which preserves the constant absolute value c of the velocity r. In 
fact, if g <^ r = 1, \q — r \ ^ 1 and the equation is just the harmonic motion f -|- r ~ q, 
of point r around q. The factor {1 — q ■ r)/{q — r)^ prevents that when we take the 
boundary value r(0) = 1, the solution does not modify this absolute value of the velocity 
of the charge. 

In the case of interaction this second equation remains the same because it 
corresponds to the definition of the centre-of-mass position which is unchanged by the 
interaction, because it only involves the U and W functions. The first equation for 
particle a is going to be replaced by dpa/dt = Fa where Pa is the corresponding linear 
momentum of each particle expressed as usual in terms of the centre-of-mass velocity 

Pa = l{qa)mqa, j{qa) = (1 - 9a)"^^^ 
and the force Fa is computed from the interaction Lagrangian ([3]) 
dLj d f dLj^ 



F 



a 



dva dt \ du„ 
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For particle 1 it takes the form: 



j., = _<,_Ii:iIi_VT^rii^ + 4f_ I (4) 



In - ■^2^ dt \2\ri - ralVl - ""i ■ *t2, 

where it contains velocity terms which behave like and acceleration terms which 
go as 1/r in terms of the separation of the charges r = jri — r2|. In this new notation 

Then the system of second order differential equations to be solved are 

Qa = {Fa - qa{Fa " Qa)) (5) 

^« = 7 — ^ \l iQa-ra), a = 1,2 (6) 

[Qa fa) 

where a is the fine structure constant once all the variables are taken dimensionless. 
For that, we take the space scale factor R = h/2mc and the time scale as T = h/2mc^. 
All terms of equation ([5]) which depend on the acceleration of the charges have to be 
replaced by the expressions of ([6]). 

It would be desirable to find analytical solutions of the above equations (JME]). 
Nevertheless we have not succeded in finding such a goal. However we shall analyse 
different solutions obtained by numerical integration. We are going to use the computer 
program Dynamics Solver [6]. The quality of the numerical results is tested by using 
the different integration schemes this program allows, ranging from the very stable 
embedded Runge-Kutta code of eighth order, due to Dormand and Prince, to very 
fast extrapolation routines. All codes have adaptive step size control and we check 
that smaller tolerances do not change the results. Another advantage is that it can be 
prepared to analyse solutions corresponding to a wide range of boundary conditions, 
automatically. 

See in figure [2] the scattering of two equal charged particles with parallel spins. 
The centre-of-mass motion of each particle is depicted with an arrow. If the two 
particles do not approach each other too much these trajectories correspond basically to 
the trajectories of two spinless point particles interacting through an instantaneous 
Coulomb force. By too much we mean that their relative separation between the 
corresponding centres of mass is always much greater than Compton's wavelength. This 
can be understood because of the above discussion about the Coulomb behaviour of 
the averaged interaction Lagrangian, if the average position of each charge is far from 
the other. For high energy interaction the two particles approach each other below 
that separation and therefore the average analysis no longer works because the charges 
approach each other to very small distances where the interaction term and the exact 
position of both charges, becomes important. In this case new phenomena appear. We 
can have, for instance, a forward scattering like the one depicted on figure [3|, which is 
not described in the classical spinless case, or even the formation of bound pairs for 
particles of the same charge, which we shall analyse in what follows. 

In figure H] we represent an initial situation for two equal charged particles with 
parallel spins such that the corresponding centres of mass are separated by a distance 
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Figure 2. The trajectories of the centres of mass and charge of two spinning 
particles with an initial centre-of-mass velocity (ja = 0.1 and a small impact 
parameter. 



Figure 3. Forward scattering of two spinning particles of the same charge with 
an initial separation 2qa{0) = 10, centre-of-mass velocity |(?a(0)| = 0.18 and a 
very small impact parameter. The two centres of mass cross very close to each 
other, with a small deviation. 
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below Compton's wavelength. Remember that the radius of the internal motion is half 
Compton's wavelength. We locate the charge labels at the corresponding points 
Va and the corresponding mass labels to the respective centre-of-mass Qa- We 
depict in part (a) the situation when the two particles have the same phase f3i = /?2- 
The forces Fa, on each particle a = 1,2, are computed in terms of the positions, 
velocities and accelerations of both charges, according to (jl]), and are also depicted 
on the corresponding centres of mass as a consequence of the structure of the equations 
(j5]). We see that a repulsive force between the charges produces also a repulsive force 
between the centres of mass in this situation. However, in part (b) both charges have 
opposite phases Pi = —(32, and now the repulsive force between the charges implies an 
attractive force between the corresponding centres of mass. If the initial situation is 
such that the centres-of-mass separation are greater than Compton's wavelength, the 
force is always repulsive irrespective of the internal phases of the particles. 



\ \ 
\ \ 
\ \ 



I I 



i ^ ^« f* I 

\ \ II 
\ \ > ' 



(a) 



(b) 



Figure 4. Boundary values for two Dirac particles with parallel spins and with 
a separation between the centres of mass below Compton's wavelength. The 
dotted lines represent the previsible clockwise motion of each charge. In (a) 
both particles have the same phase and the repulsive force between charges 
produces a repulsive force between their centres of mass, while in (b), with 
opposite phases, the force between the centres of mass is atractive. 



In figure we have another situation of opposite phases and where the initial 
separation between the centres of mass is larger but still smaller than Compton's 
wavelength. 

To analyse this situation, which is going to produce bound motions, we proceed 
as follows: We start the numerical integration by imposing the boundary condition 
that both centres of mass are at rest and located at the origin of the reference frame 
9a(0) = qa{0) = 0. For particle 2 we take the initial phase /92(0) = and for /3i we 
start with Pi{0) = and, will be increased step by step in one degree in the automatic 
process, up to reach the whole range of 27i radians. The boundary values of the variables 
ra(0) and ra(0), with the constraint |ra(0)| = 1, are taken as the corresponding values 
compatible with these phases. The whole system is analysed in its centre-of-mass frame. 
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Figure 5. (a) Another situation of two charges with opposite phases which 
produce an atractive force between the centres of mass provided they are 
separated below Compton's wavelength. In part (b), after half a cycle of the 
motion of the charges, the force becomes repulsive between the centres of mass, 
but its intensity is much smaller than the atractive force in (a) so that the 
resulting motion is also a bound motion. 



so that for subsequent boundary values these variables are restricted to qi{0) — —^'2(0) 
and qi{0) = — ^2(0). The automatic integration is performed in such a way that when 
the two particles separate, i.e., when their centre-of-mass separation is above Compton's 
wavelength, the integration stops and starts again with a new boundary value of the 
phase /9i(0) of one degree more, and the new values of the variables ra{0) and r'a(O). 
If the two particles do not separate we wait until the integration time corresponds to 
10^ turns of the charges around their corresponding centre-of-mass, stop the process, 
keep record of the phases and initial velocities, and start again with new boundary 
values. This corresponds, in the case of electrons, to a bound state leaving during a 
time greater than 10~^^ seconds. For some particular boundary values, with opposite 
phases, we have left the program working during a whole week and the bound state 
prevails. This represents a time of life of the bound state greater than 10^^ seconds. 
Leaving the computation program running for a year will only increase this lower bound 
in two orders of magnitude. The general feeling is that the bound states are sufficiently 
stable, because even the possible numerical integration errors do not destroy the stabihty. 
This process is repeated again and again by changing slightly the initial values of the 
centre-of-mass variables qa{0) and qa(0), in steps of 0.0001 in these dimensionless units 
and with /32(0) = 0, and the same procedure with /3i(0), as above. To test the acuracy 
of the integration method, we check every 10^ integration steps that the velocities of the 
charges of both particles remain of absolute value 1, within a numerical error smaller 



The whole process is repeated by changing the initial /92(0) phase to any other 
arbitrary value. We are interested to see whether different results are produced 
depending on the values of the phase difference /52(0) — /3i(0) and of the centre-of- 
mass variables qa{0) and qa{0)- We collect all data which produce bound motions, and 
find the following results: 

(i) The initial velocity of their centres of mass must be |^a(0)| < 0.01c. Otherwise the 
bound motion is not stable and the two particles, after a few turns, go off. 



than 10" 
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(ii) For each velocity |ga(0)| < 0.01c there is a range A of the pase /3i(0) = /32(0) + 7r±A 
for which the bound motion is stable. The greater the centre-of-mass velocity of 
each particle the narrower this range, so that the bound motion is more likely when 
the phases are opposite to each other. 

(iii) We have found bound motions for an initial separation between the centres of mass 
up to 0.8 times Compton's wavelength, like the one depicted in figure [5l provided 
the above phases and velocities are kept within the mentioned ranges. 

In figure [6] we show the bound motion of both particles when their centres of mass 
are initially separated qi^ = —q2x = 0.2 x Compton's wavelength, qi^ = —q2x = 0.008 
and qiy = —q2y = 0.001, (32 = and Pi = n. Now the force between the charges is 
repulsive but nevertheless, if the internal phases Pi and P2 are opposite to each other, it 
becomes an atractive force between their centres of mass in accordance to the mechanism 
shown in figure H] (b). 

This possibility of formation of low energy metastable bound pairs of particles of the 
same charge is not peculiar of this interaction Lagrangian. By using the electromagnetic 
interaction or even the instantaneous Coulomb interaction between the charges of two 
spinning Dirac particles we found in [5] also this behaviour. This bound motion is not 
destroyed by external electric fields and also by an external magnetic field along the spin 
direction. Nevertheless, a transversal magnetic field destroys this bound pair system. 




Figure 6. Bound motion of the centres of mass and charge of two spinning 
particles with parallel spins and with a centre-of-mass velocity v ~ 0.0082, for 
an initial separation between the centres of mass of 0.2 x Compton's wavelength. 

When we make the average of the position it becomes the centre-of-mass Qa 
and the repulsive force between the charges is also a repulsive force between the 
corresponding centres of mass and therefore when we suppress the zitterbewegung spin 
content of the particles there is no possibility of formation of bound pairs. 



Interaction Lagmngian for Dirac particles 



15 



Although this resuh produces a classical mechanism for the formation of a spin 1 
bound system from two equal charged fermions we must be careful about its conclusions. 
First, it is a classical description and although the range of energies which produce 
this phenomenon is a wide one it does not mean that two electrons can reach that 
binding energy. This Dirac particle is a system of seven degrees of freedom: 3 represent 
the position r, another 3 the orientation a and finally the phase /?. If we accept 
the equipartition theorem for the energy, then for the maximum kinetic energy which 
produces a bound motion mv'^/2 = 7kT/2, where k, is Boltzmann's constant and 
V = 0.01c the maximum velocity of the center of mass of each particle, then it means 
that a gas of polarized electrons (like the conducting electrons in a quantum Hall effect) 
could form bound states up to a temperature below T = 8.47 x lO^K, which is a very 
high temperature. In a second place, matter at this level behaves according to quantum 
mechanical rules and therefore we must solve the corresponding quantum mechanical 
bound state to establish the proper energies and angular momenta at which these bound 
states would be stationary. This problem has not been solved yet, but the existence of 
this classical possibility of formation of bound pairs justifies an effort in this direction. If 
the phases of the two particles are the same (or almost the same) there is no possibility 
of formation of a bound state. The two fermions of the bound state have the same 
spin and energy. They differ that their phases and linear momenta are opposite to each 
other. Is this difference in the phase a way to overcome at the classical level, the Pauli 
exclusion principle? 
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